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Deeply supercooled water exhibits complex dynamics with large density fluctuations, ice coars¬ 
ening and characteristic time scales extending from picoseconds to milliseconds. Here, we 
discuss implications of these time scales as they pertain to two-phase coexistence and to 
molecular simulations of supercooled water. Specifically, we argue that it is possible to dis¬ 
count liquid-liquid criticality because the time scales imply that correlation lengths for such 
behavior would be bounded by no more than a few nanometers. Similarly, it is possible to 
discount two-liquid coexistence because the time scales imply a bounded interfacial free en¬ 
ergy that cannot grow in proportion to a macroscopic surface area. From time scales alone, 
therefore, we see that coexisting domains of differing density in supercooled water can be no 
more than nano-scale transient fluctuations. 
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1. Introduction 

For more than two decades, it has been suggested that anomalous properties of 
liquid water reflect two distinct liquids and a low-temperature critical point at 
supercooled conditions |T]. Yet Binder has observed [2] that two-liquid criticality 
defined in terms of a divergent length scale is impossible at deeply supercooled 
conditions. Specifically, growing lengths coincide with growing equilibration times, 
and the time available to equilibrate can be no longer than the time it takes the 
metastable liquid to crystallize. In other words, metastability or instability implies 
an upper bound to the size of fluctuations that can relax in the liquid. For water, we 
argue, the bound seems to be no larger than 2 or 3 nm, corresponding to volumes 
containing fewer than 1000 molecules. 

This bound is fundamentally different than a cutoff imposed by the practicality 
of a finite simulation cell. Transient fluctuations on smaller length scales might 
seem interpretable in terms of something like a liquid-liquid transition, but the 
bound implies one can never reach large enough scales to know if that interpreta¬ 
tion is correct. The interpretation certainly seems unnecessary because reasonable 
molecular models known to not exhibit two-liquid behavior and general arguments 
independent of molecular forcefield do account for equilibrium anomalies of wa¬ 
ter HMD and nonequilibrium amorphous ices 00- 
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We will see that the principal governing parameter is the ratio of liquid’s 
metastable lifetime, tms, to its structural relaxation time, tr. If tms/tr is large 
enough, conceptions of liquid-liquid coexistence and criticality can be good approx¬ 
imations. Physically realizable examples include colloidal mixtures [5] and protein 
solutions [3|. 

While water is not among those examples, some recent simulations of water mod¬ 
els are interpreted as demonstrating two-liquid-like behavior HDHUj. Other works, 
such as Refs. m and P3], seem to discount the possibility. The disagreement is 
not an issue of the force field chosen to model water, as has been noted in Ref. he 
Rather, disagreement about two-liquid-like bistability has to do with equilibration 
or reversibility. Indeed, Ref. m shows that this bistability is reproduced by con¬ 
straining the distribution of crystal order parameter to that of the standard liquid 
state, and that this bistability disappears as the distribution is allowed to relax. 

This general result, independent of free energy sampling method, argues that 
simulations finding two-liquid behavior, for example HQHrg and HE have failed 
to control reversibility in the presence of emergent crystal order. Establishing re¬ 
versibility in molecular simulations of supercooled water is possible but difficult 
both because structural relaxation in the supercooled liquid is intrinsically slow (it 
is, after all, a glass former), and because the time scales required to relax crystal- 
order fluctuations can be orders of magnitude longer than those required to relax 
density fluctuations. (The Appendix illustrates this time-scale separation.) 

Section 2 gets to the root of the issue by providing estimates of time scales for 
supercooled water, specifically tr and tms- From those estimates, a bound for the 
correlation length of putative critical fluctuations in water is derived. Section 3 
then argues that interfacial area between domains of different densities in super¬ 
cooled water are also bounded. As such, at reversible conditions, these domains 
are transient and cannot exhibit a well defined surface tension or an interfacial 
free energy that grows as a fractional power of system size. The paper concludes 
Section 4 with a discussion of differences between supercooled water and systems 
where two-liquid behavior is possible. 

Before turning to those points, let us be clear that nothing in the analysis pre¬ 
sented below discounts the possibility of more than one nonequilibrium amorphous 
solid. In fact, estimates discussed in Section 2 make use of tools that can also be 
used for predicting the emergence of glass transitions and more than one amor¬ 
phous solid [SUE Rather than this irreversible polyamorphism, our focus herein is 
on the possibility of more than one liquid and criticality at reversible conditions. 


2. Time scales and length scales 
2.1. Regimes of supercooled water 

Figure [lja) shows a phase diagram for supercooled water. The onset temperature, 
T 0 , is the crossover temperature below which liquid dynamics is heterogeneous 
and relaxation times grow with decreasing T in a super-Arrhenius fashion. The 
liquid instability temperature, T s , is the temperature below which nanometer-scale 
domains of the liquid are unstable. Fluctuations are thus large for all T < T s , 
but because the time to reorganize the liquid grows rapidly as temperature is 
lowered, fluctuations observed over finite times can appear largest in the vicinity 
of T s . Not surprisingly, therefore, all estimates of a putative critical temperature 
for supercooled water are close to T s . 

Crystallization is avoided and glass is formed when the supercooled liquid is 
cooled fast enough and to a low enough temperature. The specific temperature for 
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p/kbar log t/1 ps p = 1 bar 




Figure 1. (a) Regimes of supercooled water, with T and p denoting temperature and pressure, respectively, 
and (b) the logarithm of the supercooled liquid’s low pressure lifetime, tms ? and structural relaxation time, 
tr. Equilibrium liquid-ice coexistence occurs at the melting temperature, T m , and the crossover to glass 
forming dynamics occurs in the supercooled liquid below the onset temperature, T 0 . For the range of su¬ 
percooled temperatures extending not too far below T m or T 0 , the rate-determining step for crystallization 
is nucleation. At yet lower temperatures highlighted with grey in Panel (b), T < T s , coarsening rather 
than simple nucleation becomes rate determining. (T 0 and T s are given in Refs. 0 and |16|.) The black 
and red lines in Panel (b) refer to theoretical predictions of tms and tr, respectively, i.e., Eqs. (16) and 
(2) of Ref. 1161 . The blue circle locates the experimental measurement of metastable lifetime by Sellberg 
et al. El. The temperature Th (= 235 K at p = 1 atm) is the homogeneous nucleation temperature given 
by Holten et al. na- It is the highest temperature at which ice crystallization proceeds quickly, i.e., where 
the metastable lifetime is of order 10 1=t:1 s. The diamonds are experimental results for tr as reflected in 
experimental measurements of viscosity na. 


this transformation depends upon the time scale at which the liquid is driven out 
of equilibrium |6j. Given this protocol dependence, it is impossible to illustrate the 
multitude of possible nonequilibrium phase behaviors in a simple two-dimensional 
plot. Nevertheless, one significant feature is captured with the line drawn in the 
figure between high-density amorphous (HDA) and low-density amorphous (LDA) 
regions. 

This line marks the p at which the T to reach the nonequilibrium time scale 
is minimum [7]. At low enough temperatures, this line relates to a first-order¬ 
like nonequilibrium transition between HDA and LDA phases. Observations of the 
transition show a large range of hysteresis with the average of the forward and 
backward transition pressures being close to that line m • In the reversible melt 
near T = T s , transient mesoscopic domains will necessarily appear either as the 
first steps in ice coarsening or as precursors to those nonequilibrium amorphous 
solids. In a molecular simulation of a small enough system carried out over a small 
enough time, those transient domains can be confused with two distinct liquids. 

For such simulations, locations for the corresponding states illustrated in Fig. 0 a ) 
differ from one model to another [163. For example, with the ST2 model of water 
used in Ref. m, the corresponding states are shifted to higher temperatures from 
those of real water by 10 to 15%. 


2.2. Metastable lifetime and structural relaxation time 

Figure 0b) shows the temperature variation of tms and tr. The lines are the 
theoretical predictions of Ref. m- At conditions where the liquid persists long 
enough to make the measurements US], the theoretical prediction for tr agrees well 
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with the experimental measurements. At lower temperatures, T < Th ~ 235K (i.e., 
in so-called “no-man’s land”), agreement with simulations at similar corresponding 
states gives further confidence in the theory that predicts the red line m- 

The theoretical prediction of tms agrees with the one experimental measurement 
of that time for T s < T < Tj-,. The metastable lifetime is a non-monotonic function 
of temperature because above T s , nucleation is rate determining in crystal forma¬ 
tion, while below T s , critical nuclei are small and plentiful, and coarsening is rate 
determining. Grey shading in Fig. [ljb) highlights the distinction between the two 
regimes. The latter regime, with its large fluctuations and coarsening, is the regime 
of interest when considering possible two-liquid behavior. 

The theoretical curve, essentially an interpolation connecting the two regimes 
m, predicts that the shortest lifetime for supercooled water is no longer than 
10~ 4 s. Experimentally, making amorphous ice by cooling requires cooling rates at 
least as fast as 10 6 K/s [2T]. Based on this cooling rate, dimensional analysis also 
gives 10 -4 s as the shortest lifetime, in harmony with the theoretical prediction. 


2.3. Largest length scale of putative critical fluctuations 

The structural relaxation time of the liquid, tr, is the time to equilibrate the liquid 
on length scale a, where a ~ 0.2 or 0.3 nm is the characteristic microscopic length 
of the liquid. Near presumed criticality, the time to equilibrate over a larger length 
scale, £, would be of order tc = tr^/cx) 2 , where z « 3 [22] . But as Binder notes 
[2], < tms because the liquid will not resist crystallization for times longer than 

r M s- Accordingly, £/a < (ths/Vr,) 1 / 3 - 

This bound giving the largest correlation length for critical-like fluctuations is 
applicable in the regime where criticality is imagined to occur, T < T s . Figure [I] 
shows that in this regime at ambient pressures, tr and tms grow with decreasing 
temperature, and the ratio is tms/tr 1$ 10 3 throughout. Based upon the theoretical 
predictions, we expect tms/'Tr 1$ 10 3 to remain true at the slightly elevated pres¬ 
sures that are sometimes identified with the putative liquid-liquid critical point, 
e.g., p ss 0.5kbar [TH] , 

Specifically, according to the theory [16], increasing pressure up to 2 kbar will 
decrease the ratio of tms/tr, but by no more than 10 to 20 %. The decrease is due to 
the extent by which increasing pressure increases (<9 tr JdT) p [7] and decreases the 
enthalpy of fusion [23]. The decrease in enthalpy of fusion causes tms to decrease 
because the solid-liquid surface tension is proportional to the enthalpy of fusion 
[24] . a proportionality that has been tested for a molecular model of water [25] . 
Experimental tests have not yet been performed to test the prediction of tmsAr 
at elevated pressures. 

Accepting the prediction, as the theory compares favorably with experiment at 
ambient pressures, it follows that by cooling the system down to T s , the corre¬ 
lation length of density fluctuations (which may have the physical interpretation 
of precursor effects of the LDA-HDA nonequilibrium transition at still lower tem¬ 
peratures) can increase by at most one order of magnitude, i.e. £ < 2 or 3 nm. 
Phenomena with such a small largest length scale would seem to be poorly approx¬ 
imated by criticality, which is defined by a diverging correlation length, although 
in a computer simulation context such a modest increase of the correlation length 
might be easily mistaken as a signature of a true transition. 
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3. Interfaces between finite domains in the metastable liquid 

With a similar argument, one can also conclude that different liquid domains at 
metastable conditions are bounded in size. To do so, bear in mind that where 
nucleation is rate determining, the liquid’s lifetime decreases with increasing system 
size (i.e., the bigger the system the more opportunities for a critical nucleus to 
appear). In contrast, in deeply supercooled conditions, where coarsening is rate 
determining, the metastable lifetime increases with increasing system size (i.e., the 
bigger the system the more likely to have misaligned crystal domains). Dependence 
on system size in the latter case is weak, growing as a fractional power [26] . 

The range of temperatures considered in Figure [T](b) , T < T^, is the crossover 
between the two regimes. Given the opposing system-size dependence, we can ex¬ 
pect the system size dependence to be especially small in the crossover regime. 
Whether negligible or weak, the implication is that there can be no surface tension 
for interfaces separating domains of different metastable states. This implication 
follows from log(rMs) providing an upper bound to any free energy barriers that 
might separate distinguishable states in the metastable liquid. In particular, two 
coexisting domains separated by a stable interface will interconvert on a time scale 
ri n t, and this time is proportional to exp(/3A F), where A F is the interfacial free 
energy and (5 is reciprocal temperature. To be observable, this time must be smaller 
than tmsQ 

As a result, the free energy barrier for interconverting coexisting domains of 
low-density liquid and high-density liquid cannot be an interfacial free energy that 
would grow (unbounded) as a fractional power of system size. Rather, this interfa¬ 
cial free energy grows at most logarithmically with system size, which implies that 
different forms of supercooled liquid water can be at most transient and finite in 
size. 

To assign numbers to this bound, we can use transition state theory for a diffu¬ 
sive barrier, i.e., Ti nt ~ Aexp(/3AF), where A ss Tc^tr/to) 1-61 . Here, r 0 ~ 1 ps is 
the characteristic microscopic time |29], and 9 ~ 1/4 is the decoupling exponent re¬ 
lating reorganization time to diffusion ES EH. Thus, because Ti nt < tmSj we arrive 
at the bound f3AF < (1/4) ln(rR,/r 0 ) + R^ths/tr)- Applying this formula in the 
regime where different liquid domains appear, i.e., T < T s , we arrive at /3AF < 10. 
This microscopic bound to A F implies that away from criticality, interfaces of area 
larger than microscopic scales cannot be stable. Indeed, recent simulation data is 
consistent with this bound, with free energy estimates exhibiting a weak system 
size dependence before saturating at a value of j3AF less than 10 [TP] , 


4. Discussion 

The arguments we present herein rely on known time scales for relaxation in super¬ 
cooled water and established theory relating time scales to equilibrium correlation 
lengths and free energy barriers. While we focus on water, the arguments apply 
equally well to other supercooled liquids. For some such materials, it is possible 
for tms/tr. to be sufficiently large that criticality and two-liquid coexistence can 
be good approximations. As noted in the Introduction, an assortment of simula¬ 
tion models and colloidal suspensions behave in this way. See, for example, Refs. 

ME2E31. 


1 This argument presumes that reversible two-liquid behavior requires interconversion without involvement 
of the crystal phase. Were ice to be an intermediate for interconversion between a low density liquid and 
a high density liquid, the bound discussed herein would not apply. 
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log 10 t/MCS log 10 (/MCS 

Figure 2. Relaxation functions for a variant of the ST2 model of water at p = 2.2kbar and T = 235 K as 
a function of Monte Carlo steps (MCS). Red and blue lines are, respectively, the p and Qq autocorrelation 
functions in two different different windows sampled during free energy calculations of Ref. [13] . The 
average value of p and Qq in a specific window is p and Qq. Random oscillations about zero at the largest 
times are the results of autocorrelating over finite times, typically between 50 to 100 times that for the 
Q6-correlation function to reach 0.1 of its initial value. Notice two (or more) step relaxation for p, and 
that for small Qq, the long-time relaxation times of p increase and approach those of Qq as p decreases. 


While divergent scales and two-liquid coexistence are impossible for metastable 
liquids in principle, nothing prohibits such behavior at conditions where a fluid 
is stable with respect to the crystal. Recent work by Smallenburg et al. [33] is a 
case in point. That work examines a class of models with nominally tetrahedral 
patchy particles. Stable liquid-liquid phase transitions do occur for some members 
of that class, those with parameters that suppress tetrahedral ground states in 
favor of BCC orderings. When parameters are changed to favor water-like models 
(i.e., locally tetrahedrally ordered), liquid-liquid behavior becomes metastable and 
therefore not observable beyond microscopic scales. 
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Appendix 

Time-scale separation 

Separation between time scales for fluctuations in density, p, and time scales fluc¬ 
tuations in global crystal order, Qq, are self evident because density fluctuations 
can be observed over periods of time where liquid water persists at supercooled 
conditions. The nature of these differing scales can be examined with simulation 
by studying autocorrelation functions for these variables. Doing so, however, is 
not straightforward because distributions of supercooled-liquid micro states are 
not stationary. On the other hand, distributions are stationary for each of the 
constrained ensembles sampled during during the Monte Carlo calculations of the 
reversible-work surface, F{p,Qo). 
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While it is difficult to assign an absolute physical time to the steps of a Monte 
Carlo trajectory, one may nevertheless examine relative relaxation times with the 
autocorrelation functions obtained within a given sampling window, 6, as a function 
of time t in units of Monte Carlo steps. Each different value of b coincides with a 
different window and thus a different value for the average p and Qq within that 
window. 

Figure [2] shows such relaxation functions obtained in this way for p and Qq 
obtained from data assembled for Ref. [13] for the ST2 model, 

(8p(0)6p(t)) b {8Qa(0)8Q 6 (t)) b 

{{5 P f) b an m,f) b ’ 

respectively. Here, (• • • ) b denote ensemble average with the biasing potential used to 
confine configurations to the 6th window of p-Q§ space in a free energy calculation. 
The fluctuations, bp and 5Qq, are deviations from their respective means in the 
6th window. 

The correlation functions for the two windows considered in Fig. [2] are repre¬ 
sentative of what we find over a swath of more than 100 separate constrained 
ensembles covering a broad variety of relaxation behaviors. The time-scale sepa¬ 
ration between p and Qq is especially clear at the higher density characteristic of 
early-stage coarsening. The time-scale separation diminishes at the longer times 
but remains significant at the shorter times for the lower density characteristic of 
ice. 
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